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To clarify the evolutionary relationships among 
betavoronaviruses that infect bats, we analyzed samples 
collected during 2010-2011 from 14 insectivorous bat 
species in China. We identified complete genomes of 2 novel 
betacoronaviruses in Rhinolophus pusillus and Chaerephon 
plicata bats, which showed close genetic relationships with 
severe acute respiratory syndrome coronaviruses. 


he 2003 outbreak of severe acute respiratory syndrome 

(SARS) was caused by a novel betacoronavirus and 
rapidly spread globally, causing ~8,000 cases and nearly 
900 deaths (/,2). In June 2012, a novel betacoronavirus 
(called human coronavirus EMC [HCoV-EMC)]) also was 
isolated from the sputum of a patient from Saudi Arabia 
who died of pneumonia and renal failure (3). Similar vi- 
ruses were detected in 2 additional patients who had severe 
pneumonia in Qatar in September 2012 and in Saudi Arabia 
in November 2012 (4,5). The clinical picture was remark- 
ably similar to that of SARS and illustrates the epidemic 
potential of a novel coronavirus (CoV) to threaten global 
health. SARS-CoVs and HCoV-EMC were suspected of 
spreading from bats to humans because these CoVs were 
most closely related to bat CoVs (/,4). To clarify the evo- 
lutionary relationships among betavoronaviruses that in- 
fect bats, we analyzed samples collected during 2010-2011 
from 14 insectivorous bat species common in 8 provinces 
in China. 


The Study 

We obtained pharyngeal and anal swab specimens 
of 414 insectivorous bats. Samples of each species were 
pooled and then processed with a viral particle—protected 
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nucleic acid purification method (6). The extracted RNA 
and DNA were amplified by sequence-independent 
PCR. The amplified viral nucleic acid libraries of the bat 
species were then sequenced with the Illumina/Solexa 
GAIT sequencer (Illumina, San Diego, CA, USA). Those 
reads generated by the Illumina/Solexa GAIT] with length 
of 80 bases were directly aligned to the protein sequences 
in the National Center for Biotechnology Information 
nonredundant protein database by the blastx program in 
the BLAST software package, version 2.2.22 (www.ncbi. 
nlm.nih.gov/blast) with parameters “-e le-5 -F T -b 10 
-v 10.” No assembly was performed before alignment. 
Sequence similarity—based taxonomic assignments were 
conducted as described (7). We found 1,075 reads of 
betacoronavirus in Rhinolophus pusillus bats in Shaanxi 
and 92 reads of betacoronavirus in Chaerephon plicata 
bats in Yunnan. 

We estimated the approximate locations of those 
reads on the CoV genome and their relative distances on 
the basis of alignment results exported with MEGAN 4— 
MetaGenome Analyzer  (http://ab.inf.uni-tuebingen.de/ 
software/megan/). The located reads were then used for 
reads-based nested PCR to identify genomic sequences. 
We established the complete genome sequences of 2 
betacoronaviruses (Bat  Rp-coronavirus/Shaanxi2011 
and Bat Cp-coronavirus/Yunnan2011), which are 29,484 
nt and 29,452 nt, respectively. The G+C content of Bat 
Rp-coronavirus/Shaanxi12011 and Bat Cp-coronavirus/ 
Yunnan2011 is 41.6% and 40.9%, respectively. 

We conducted complete genome comparison and 
phylogenetic analysis on the basis of polymerase and 
spike protein. Pairwise genome sequence alignment 
was conducted by using EMBOSS Needle software 
(www.ebi.ac.uk/Tools/psa/emboss_needle/) with default 
parameters. The overall nucleotide sequences between 
Bat Rp-coronavirus/Shaanxi2011 and Bat Cp-coronavirus/ 
Yunnan2011 indicated 88.7% nt identity. They shared 
87.4%-89.5% nt identity with SARS-CoV, 88%-89.9% 
nt identity with the bat SARS-like CoV (bat SARS- 
CoV Rm1), and 87.6%-89.6% nt identity with the civet 
SARS-like CoV (civet SARS-CoV SZ16). On the other 
hand, comparison between the betacoronavirus genomes 
and human betacoronavirus (HCoV-OC43) showed 
only 49.9%-50.4% nt overall identity, whereas the 
betacoronavirus genomes and HCoV-EMC showed 52.1% 
nt overall identity. 

The RNA-dependent RNA _ polymerase (RdRp, 
the 12th nonstructural protein codified to open reading 
frame la,b) is a highly conserved gene of CoVs, which 
is frequently used for phylogenetic comparison (8,9). 
MEGAS.0 (www.megasoftware.net) was used to construct 
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DISPATCHES 


the phylogenetic trees on the basis of the nucleotide 
sequences and deduced amino acid sequences. First, we 
used the MUSCLE package and default parameters (www. 
megasoftware.net/) to construct the alignment. The best 
substitution model was then evaluated with the Model 
Selection package implemented in MEGAS. Finally, we 
used the maximum-likelihood method with an appropriate 
model to process the phylogenetic analysis with 1,000 
bootstrap replicates. We constructed a phylogenetic tree 
based on the nucleotide sequences of the RdRp gene 
to show the evolutionary relationship between these 2 
betacoronaviruses and other CoVs (Figure 1). Reference 
CoV genome sequences were downloaded from GenBank 
and aligned with the fragments of the newly discovered 
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Figure 1. Phylogenetic tree of novel betacoronaviruses based 
on the nucleotide sequence of the RdRp gene. The following 
coronaviruses (CoVs) and GenBank accession numbers were used: 
bat severe acute respiratory syndrome CoV Rm1 (bat SARS-CoV 
Rm1; DQ412043), bat SARS-CoV Rp3 (DQ071615), bat SARS- 
CoV Rf1 (DQ412042), bat SARS-CoV HKU3 (DQ022305),SARS- 
CoV isolate Tor2/FP1—10895 (SARS-CoV Tor2; JX163925), SARS- 
CoV BJ182-12 (SARS-CoV BJ182; EU371564), SARS-CoV 
(NC004718), civet SARS-CoV SZ3 (AY304486), civet SARS-CoV 
SZ16 (AY304488), bat CoV HKU9 (BtCoV-HKU9; EF065513), bat 
CoV HKU4 (BtCoV-HKU4; EF065505), bat CoV HKU5 (BtCoV- 
HKU5; EF065509), human betacoronvirus 2c EMC/2012 (HCoV- 
EMC; JX869059), human CoV OC43 (HCoV-OC43; NC005147), 
HCoV-HKU1 (NC006577), bat coronavirus HKU2 (BtCoV-HKU2; 
NC009988), bat coronavirus 1A (BtCoV-1A; NC010437), HCoV- 
229E (NC002645), HCoV-NL63 (NC005831), bat CoV HKU8 
(BtCoV-HKU8; NC010438), scotophilus bat CoV 512 (BtCoV-512; 
NC009657), avian infectious bronchitis virus (IBV; NC001451), 
beluga whale CoV SW1 (BWCoV; NC010646). Scale bar indicates 
genetic distance estimated by using TN93+G+I model implemented 
in MEGA5 (www.megasoftware.net). 
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Shaanxi2011 and Bat Cp-coronavirus/Yunnan2011 were 
highly similar, sharing 93.1% nt identity. The phylogenetic 
analysis demonstrated that betacoronaviruses and the bat 
SARS-like CoVs in our study are clustered (93.1%—93.4% 
nt identity) and are close in distance to SARS-CoVs 
(92.9%-94.8% nt identity) and civet SARS-like CoVs 
(93.1%—-94.8% nt identity) but that bat CoV (BtCoV- 
HKU9) and HCoV-OC43 are placed among the relatively 
distant groups (65.8%-65.9% and 62.9%-63.5% nt 
identities with the betacoronaviruses, respectively). 
Therefore, collectively we called these betacoronaviruses 
and bat SARS-like CoVs the bat SARS-like cluster of 
CoVs. Bat Rp-coronavirus/Shaanxi2011 and Bat Cp- 
coronavirus/Yunnan2011 showed little genetic similarity 
(<66.2%-—67.3% nt identity) to HCoV-EMC. 

The spike proteins of CoVs are responsible for receptor 
binding and host species adaptation, and their genes 
therefore constitute one of the most variable regions within 
CoV genomes (/0,//). The phylogenetic tree based on the 
amino acid sequences of spike protein (Figure 2) suggests 
that the selected betacoronaviruses were mainly divided 
into 5 clusters: SARS cluster; bat SARS-like cluster; civet 
SARS-like cluster; human betacoronavirus cluster; and 
EMC cluster. Bat Rp-coronavirus/Shaanxi2011 and Bat 
Cp-coronavirus/Yunnan2011 shared 89.4% aa identity 
in spike proteins, which consisted of 1,240 aa and 1,241 
aa, respectively. The spike proteins of the CoVs in our 
analysis have 89.8%-92.7% aa identity with those of 
bat SARS-like CoVs, with substantial similarity in the 
receptor-binding domain. The close relationship also was 
observed with the SARS-CoVs (79.2%-—79.4% aa identity) 
and civet SARS-like CoVs (78.9%—79.1% aa identity). In 
contrast, the human betacoronaviruses and EMC cluster 
formed separate clusters distinct from SARS-related CoVs 
that showed only 27.8%-29.4% aa and 28.8%-30.5% 
aa identities with the betacoronaviruses, respectively, 
in our analysis. The genome sequences reported here 
have been deposited into GenBank (accession nos. 
JX993987-JX993988). 


Conclusions 

The recent fatal human infection caused by HCoV- 
EMC has boosted interest in the discovery of novel 
CoVs in humans and animals. HCoV-EMC is a novel 
betacoronavirus, and its closest known relatives are 
BtCoVs HKU4, and HKUS, which have been detected in 
Hong Kong only in bats (/2), the same animal from which 
SARS is believed to have originated. Bats are increasingly 
recognized as natural reservoirs of CoVs and may serve as 
intermediate hosts for interspecies transmission of SARS- 
CoVs (/0,/3). Different bat populations from various 
countries harbor diverse CoVs that have a high frequency 
of recombination and mutation rates that enable them to 
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Figure 2. Phylogenetic tree of novel betacoronaviruses based on the 
deduced amino acid sequence of spike protein. SARS, severe acute 
respiratory syndrome; CoV, coronavirus; HCoV, human CoV; BtCoV, 
bat CoV; BWCoV, beluga whale CoV; IBV, avian infectious bronchitis. 
Scale bar indicates genetic distance estimated by using WAG+G+I+F 
model implemented in MEGA5 (www.megasoftware.net). 


adapt to new hosts and ecologic niches (/4,/5). Therefore, 
continuous studies of CoVs from different bat species and 
different countries would help better prevent the new 
global pandemics resulting from novel viral infection. 

We detected and _ characterized 2 _ novel 
betacoronaviruses—Bat Rp-coronavirus/Shaanxi2011 in 
R. pusillus bats and Bat Cp-coronavirus/Yunnan2011 
in C. plicata bats—in China. The high similarity shown 
by phylogenetic analysis confirmed the close genetic 
relationship among the CoVs (SARS-like CoVs and 
SARS-CoVs) that we analyzed. In contrast, Bat Rp- 
coronavirus/Shaanxi2011 and Bat Cp-coronavirus/ 
Yunnan2011 showed little genetic similarity with human 
betacoronaviruses and HCoV-EMC. Although several 
CoVs are found in horseshoe bats (Rhinolophus spp.), to 
our knowledge, the SARS-like CoVs in R. pusillus and 
C. plicata bats in China have not been identified. The 
description presented here will further the understanding 
of CoVs distribution in different bat species found in 
human habitats and provide clues for rapid response to 
potential public health threats. 
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